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For Gagen-Poiseuille flow, we show that exponential instability (to extremely small, axially 
symmetric disturbances represented by Galerkin's approximation) is possible only if there exists 
conditionally periodic variability of the disturbances along the pipe axis when the threshold 
Reynolds number Re.r(p) depends on the ratio pof two longitudinal periods. Absolute 

minimum Re ( ,(p) « 448 (for p « 1.53..) is obtained that corresponds to the observed conditions 

of transition from the laminar resistance law to the turbulent one and Tollmien-Schlichting 
waves exciting in the boundary layer. 
PACS: 47.20.Ft, 47.27.Cn, 47.27.nf 

1. Fundamental and applied problem of understanding of the mechanism of turbulence 
development for the Gagen-Poiseuille (GP) l flow remains mysterious because of the paradox of 
the linear stability of this flow to extremely small magnitude disturbances for any large Reynolds 
number values Re [1-4]. Up to now, to bypass this obvious paradoxical contradiction to 
experiment, a consensus is made to allow only non-linear mechanism of the GP flow instability 
to a disturbance with a sufficiently large finite magnitude [5-8]. This assumption is usual (see [3, 
4]) based on the very special interpretation of the experiments in which the substantial increase 
of the threshold Reynolds number Re.- value (defining transition from the laminar to turbulent 

flow state) is achieved by increasing the smoothness of the streamlined pipe surface. In such an 
interpretation, only correlation between smoothness increasing and resulting decreasing of the 
initial disturbance magnitude is taken into account. Meanwhile, noted already by O. Reynolds 
[1] extremely high sensitivity of the value of Re>r to the initial disturbance does not exclude 

possibility of affecting on Re., not only magnitude but frequency characteristics of the 

disturbances (resulting from non-ideal smoothness of the streamlined surface) as well. Actually, 
for example, in the experiment [9], it was found that under the fixed magnitude of artificially 
excited disturbances, GP flow instability develops only in the narrow disturbances frequency 
range. 

Herein, we show that possibility of the linear absolute instability of the GP flow in the general 
case also is defined by the value of the additional to Re control parameter p , which 
characterizes frequency properties of the disturbances and affects on the value of Re., {p) 

independently from the magnitude of the initial disturbances. Such a parameter p introducing is 
performed below on the base of pointed by O. Reynolds [1] (and by W. Heisenberg, see [4, 6]) 



1 GP flow is a laminar, stationary flow of the uniform viscous fluid along the static, direct and the length unbounded 
pipe with the round the same along the pipe axis crosscut. 

2 For special GP flow modifications (in the pipe with elliptic crosscut [10], or in the rotating pipe [11], or in the case 
of the flow with the finite size particles transported by it[12]), there initially exists natural additional to Re control 
parameter, and the paradox of linear stability is absent. 
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dissipative GP flow instability mechanism related with the action of molecular viscosity v in 
the proximity of the solid boundary. According to [1], the mechanism manifests itself as a 
spontaneous one-stage appearance for Re>Re.r of the vortexes having character size / , 

«.. which does not grow any more contrary to expectations with the growth of the velocity 
magnitude [1]». As a result, in addition to the character disturbance scale related to the pipe 
diameter, 2R , due to near-boundary vortex-generating viscosity action, there exists a new 
additional scale / ,which is proportional to the size of some part (which is really activated only 

for v^O) of the pipe boundary. It, together with R , can define frequency parameters of the 
initial disturbances, e.g., the longitudinal along the pipe axis (axis z ) spatial periods. The ratio 

/ 

of the periods, p = — , as it is shown below, is a new additional to Re parameter defining the 
2R 

GP flow instability threshold to the conditionally periodic on z extremely small vortex 
disturbances. Such representation of the disturbances structure meets in observed conditionally 
periodic Tollmien-Schlichting (TS) waves development of which (caused also by near-boundary 
action of the molecular viscosity) precedes explosive emerging of the turbulence in the near- 
boundary layer [17-19]. Moreover, already in [2, 20], it was noted that usually considered in the 
linear stability theory strictly periodic on z structures of the disturbances field obviously do not 
match those observed in the experiments. We show further, that the assumption of the strict 
longitudinal (by z ) periodicity of the disturbances may result in the paradoxical conclusion on 
linear stability of the GP flow for any Re with Re.r — »oo. The same time, with the simplest 

Galerkin's approximation, we get the finite value of the absolute minimum Re,r « 448 (for 

p « 1.53..), that is the observed characteristic threshold value Re for transition from the laminar 

resistance law to the turbulent one [2, 21] and for the conditions of the TS waves exciting [19]. 

Qualitatively the proposed theory conclusions agree with the experimental data for the flows in a 

pipe [22-24]; also, we have conducted comparison of our conclusions with conclusions of the 

stability theory (Schlichting and Lin's), and with experimental data on the stability of laminar 

near-boundary layer [25]. 

2. Let's consider the well known (see [4]) representation of the GP flow in the cylindrical 

r 2 p2 dpp. 

coordinate system (z, r, <p) : V Qr = V Q(p = 0,V Qz = V WBX (1-—^), where Enax = 4 £ , fluid 

dpn 

density p — const , -^- is a constant value of the pressure Pq gradient along the axis of the 

pipe with radius R , and v is the fluid kinematic viscosity coefficient. 

In the axially symmetric case (i.e. for extremely small disturbances which are not depending on 
<p), linear instability of the GP flow may be defined only by the value of the tangential 
component of the disturbance velocity Vm , which before coming to the non-linear evolution 
stage does not depend on other disturbance fields and satisfies the following equation: 

dK sv a v m 

-f + F 0z (r)-^ = v(A^-4), (1) 

dt oz r 

where A is the Laplace operator. It follows from (1) that exponential growth of V with time, 

not depending on the reference system choice (i.e. absolute instability of the flow [2, 4]), 
evidently is not realizable as for v - so for the case of strictly periodic or localized by z 



3 Such a mechanism is realized in the systems having disturbances with the negative energy [13-16] as , for the 
threshold generation of the vortexes (rotons) in the stream of the superfluid helium in capillary [13]. 



functions Vm in (1). It is easy to verify the latter, if to multiply (1) by Vm and integrate the both 

sides of (1) by z in the limits of the period of the function Vm (or from -oo to +co in the case 

of the localized by z function Vm ), when the contribution of the second term in the left hand 

side of (1) is zero. 

Thus far, only for conditionally periodic by z (see [2]) functions and only for v^Oin (1) it is 
possible to expect realization of the absolute dissipative instability of the GP flow for above the 

'max 

threshold Reynolds numbers Re > Re .. , where Re = . 

in v 

Let's find solution of (1) in Galerkin's approximation (see [3]), describing dependence of Vm 

on r as follows 

V^V^ASz^J.U,/-) , (2) 

n=\ & 

where Vy in (2) for any A n meets necessary boundary conditions on r ( V^ < go for r — , and 
V(p =0for r = R), because J-y is the first-order Bessel function, and y-y are its zeroes, i.e. 
J-\(j\ ) = 0for any n = l,2,.. For coefficients A n , from (1), (2), we get the following system of 
equations in the dimensionless form: 

^ + JlA. - ^ + Re ±P m ^ = 0, (3) 

OT OX n -\ OX 



t V z 

where m = \,2,..J\[,T = —^,x = - B . In (3), constant coefficients, P nm , are of the form 

R 2 R 

2 ' 
P nm =- f Y^-ldyy(l-y 2 )J l (j\ n y)J 1 (j lm y) , (4) 

where J 9 is the second-order Bessel function, and the linear on y term under the integral yields 

' . For TV = 1 in (3), the last 

0,n^m 

term may be excluded by the Galilean transformation, and hence for N -I, there is no 
possibility for the absolute instability of the GP flow. In relation with this, we shall further 
consider (3) for the simplest non-trivial case, N — 2, that allows already to eliminate the GP flow 
linear stability paradox, and leads to the conclusions quantitatively agreeing with the 
experimental data (see section 4) [22, 24]. 

3. Observed in the experiment field structures do not agree with the strictly periodic along the 
pipe axis disturbances (see herein above, and [2, 20]). Moreover, in [20], it is noted that different 
radial modes (defining dependency on the radial coordinate) correspond to different variability 
periods along the pipe axis. The latter can be modeled with the use in (3) of different periods 
along the pipe axis for the modes with different index value m. Such a requirement corresponds 
to introducing for each of these modes of its own, independent from other modes, periodical 
boundary condition. Hence, it is necessary to use the Galerkin's approximation for a dependency 
of the solution of (3) on x relevant to the mentioned above boundary conditions. Let, thus far, in 
(3), for example for N-2, magnitudes A and Aj have the form of the running waves with 

different periods along the pipe axis: 

A x = A x Q e A r+aiaxx ? ^ = A ^ r+i2nf5x ? (5) 



where A, q and Ay\ are constants, and the additional to Re control parameter is defined as 

a l 

p = -n for all a and B , when, in particular, p = — = -^— . In Galerkin's approximation (for the 

solution of the system (3) depending on x in the form (5) for a ■*■ /3 generally), we get a uniform 
system for A,q and Ay^ , condition of solvability of which for the finite A>q and ^n ai l° ws to 

define the value of the exponent X = \+ iX 2 depending on dimensionless parameters Re , p , 
and B (see (A.l) and (A.2) in Appendix). 

Vortex structure of the considered in (2), (5) disturbances is related to non-zero components 
(longitudinal and radial) of the vorticity field, and zero value of the integral value of the Lamb's 
vortex momentum density ,which for the main undisturbed GP flow has longitudinal component 

pV 
only equal to — 5iM . This structure has finite value of the integral helicity for the resulting flow 

which is a linear superposition of the main GP flow (having non-zero the longitudinal 
component of the velocity field and the azimuthal component of the vorticity only) and of the 
disturbance field. 

Usage according to [20] of the form (5) to get solution of the system (3) just in the Galerkin's 
approximation (instead of usually applied exact solution of the system (3) for the radial modes 
having one and the same longitudinal period) may be also interpreted as a result of the possible 
modeling of the always existing stochastic heterogeneities of the pipe internal surface. 

Linear exponential instability condition with A. > Oin (5) has the form of (A.3). For Re »1, 

it may be reduced to more simple (A.4).Then for /? from (A.5), minimizing Re , in (A.4), it 

defines the following linear exponential instability condition for GP flow (providing a lower 
boundary for an exact value of Re , , determined from (A.3)): 

Re>Re ( ^ 2 "-^f\ (6) 

p 12 p 21 p 2 \s\ 

n \ S 
where S* = sin ;zp sin — sin n{ — h p) , B = (pP - P ) , 

p p | S | u ll 

F = (yl 2 +r l,l )il + p2)A2+(r l,2 -rl x )(\-p 2 )B 2 +2AB(yl a -p 2 /^) , 

2 2 4^^/c^ + i) 

A -B — for P. -,,P ni P.~,P~. from (4) , for arbitrary value of p 

n (\-p) 

9 
and for A > . 

2 
Only for A < 0, it is necessary to use instead of (A.4) and (6) the general condition (A.3), for 

which A > in (5). 

In (6), infinite value Re , — > oo takes place only for those finite values of p for which S = , 



1 k + \±yl(k + \) 2 -4 

i.e. for p = p 1 =k,p = p^ =—, and p = p , = for all integer 

* _ k VK+1 2 

k 
£ = 1,2,... For p, being inside the intervals of changing of p between any two consecutive 



values of p 1 ,p 1 ,p , , value of Re , in (6) is a function of p having one local minimum 

* j_ -Jk + 1 '" 

ifc 

inside each of the pointed out intervals (see Fig. l,b). So far, the absolute minimum value, 

Re^ 1 « 442 , in (6) is reached for p « 1.53.., close to the "golden" ratio p = l±^l « 1.618 .. 

(i.e. the limit of the infinite sequence of ratios of two consequent Fibonacci numbers 1, 2, 3, 5, 8, 
13, 21, ..). For the same/? , from the exact condition (A.3), we get close to it absolute minimum 

value Re , « 448 (see also Table in Appendix where conclusions based on (A.3) and (6) are 

compared). 

4. Stated value Re , « 448 corresponds to the interval of Re e 300-^500 noted in the 

experimental observation of the threshold transition from the laminar resistance law (for the flow 
in the pipe) to the turbulent one [2, 21], and of the TS waves exciting in the near-wall region of 
the boundary layer [19]. Observed in [1] and other experiments (see references in [23, 24]) 

hyper-sensitivity of the value of Re , to the initial disturbances, actually, corresponds to the 

obtained in (6) dependency of Re , on p when, for example, Re , in (6) changes 

approximately 600 times with the change of /?from 0,1186 to 0,1118 only. Neighboring local 

minima of Re , in (6) may also significantly differ each from the other, when for p « 2.23 , we 

have in (6) Re , ~ 1982, and for p « 3.86, we already have Re , ~ 84634. Enlarged fragments 

of the neutral stability curve corresponding to the condition (6) (see Fig. 1, b), are represented in 
Fig. 1, a, as dependence of the value of l/2p on Re. They are plotted on the taken from the paper 
[25] figure (see Fig. 12 in [25]) representing theoretical (Lin's and Schlichting's) neutral stability 
curves and corresponding to them experimental data related to the determining of the instability 
threshold in the boundary layer. 

Considered in (2), (5) conditionally periodic on z structures of the initial disturbance field 
V qualitatively agree with observed in [23] wave-like changes of V . The latter are especially 

clearly seen in the proximity of the turbulence decay threshold for Re ~ 1750 when the ratio of 

8 13 

visible in Fig. 5d) from [23] character periods with p « - and p « — is close to the pointed out 

5 8 

above "golden" p . In Fig. 2a), taken from [24] (see Fig. 4 b [24]), it is giving the dependency 

on Re for the constant velocity of the rear front of the observed in the flow in the pipe [24] 
turbulent spot together with the experimental data [22] related to the given in Fig. 2a) interval of 
the values of Re (dark triangles for the rear front of the turbulent spot, and white triangles for the 
leading front). Also, results following from the non-linear theory [8] and for the phase velocity 

IrJ V 

v = — L max ( scaled by the mean flow GP ve i oc ity y = _max_ ^ defined in ( A 6 ) from the 
P InfiRs, m 2 

representation of /L in (A.2) for the neutral curve (i.e. when A. 1 = ) are presented there. 

Hence, a theory based on the simplest Galerkin's approximation with TV = 2 in (3) already 
yields agreeing with experiments [22, 24]. The same time, its worth to mention that linear 
exponential instability in (6) turns out to be possible also for rational p not equal to p = p k or 

p , , but for irrational p-p , , k = 1,2,.., vise versa, Re , — > oo in (6). It means that these 

1 V& + 1 tn 



conclusions based on the Galerkin's approximation with N — 2 may in this respect differ from 
the made above (in the section 2) general conclusions on the base of the original equation (1). 
Future GP flow stability analysis for N > 2 in (3) shall clarify the reason for such a discrepancy. 

Obtained results allow to fill a gap in the non-linear theory [7, 8] in which up to now instead 
of the linear exponential instability it was necessary to consider a stage of the seed algebraic 
instability (where small initial disturbance may grow only locally in time tending to zero with 

We are grateful to S.I. Anisimov, E. A. Novikov and N.A. Inogamov for useful comments 
and interest to the work. 

Appendix 

1. From (3) and (5) for TV = 2 in Galerkin's approximation (representing disturbance 
dependency on the longitudinal coordinate z when generally a ^ /? in (5)), one gets for 

a = Aj + i/L 2 '. 



1/2 
1/2 
1/2 1/2 



\ = ~ r \,\~ 47r2 ^ 2p2 ~2 (a l ± J^ D ( } ' (A>1) 

A 2 = -2n/3pP n Re- i(a 2 ± -j=D£ ) , (A.2) 



2 2 2 

where Z) = d' +1 , D =d' -I , 1 = a, -a„ + Ac^ Re 1 

a l =rl 2 -rl l +4x 2 j3 2 (l-p 2 ), a 2 =2;&Re(P 22 -pP n ), 

2 Re2p 21 P MP 2 P 2 2 2 2 

Re 2 = 2J_i^ ^ , = l z +4(a a -2Re fd ) , d =-AS, 

1 (\-p) 2 U Z til 

c 1 = -4Sctgn(p H — ) , and S is defined in (6). 
1 p 

Condition X x > leads to the inequality 



d 



(aRe+bV >c + — - (A.3) 

Re 2 

AP P p 2 p 2 S 
where a = Z[ [Z , b = ifi{P - pP )a , 

(1-/0 

d = a 2 (y 2 l+ 4x 2 j3 2 p 2 )(y 2 2 +4x 2 j3 2 ) , a 3 =y 2 2 +y 2 l+ 4x 2 j3 2 (l + p 2 ), 

2 1 

P \1?2\ P Sct § 7r (P + -) 

c = a 2 /3 2 (x 2 (P -pP ) 2 -4 P-). 

* (1-/0 

2. In the limit of Re » 1 , inequality (A.3) for c > is reduced to the inequality 

Re>Re^(/?) = — ^1 . (A.4) 

In (A.4), function Re , (/?) has the minimal value (given in (6)) when 



P = Pc 



1 



^r^rh+Bir^-rh 



A(\ + p 2 ) + B(\-p 2 ) 



In 

where A and B are defined in the main text (see (6)) for A^ > . 
3. On the neutral curve 

velocity V p I V cp is as follows 



(A.5) 



3. On the neutral curve with X.-O (i.e. when equality takes place in (A.4)), the phase 



VJV c =pP u +P 22 ±(- 



u 2 



")i 



(A.6) 



where B = /? 12 (Re, p) corresponds to the replacement of inequality by equality in (A.4). With 
such replacing of inequality by equality in (A.4), one gets quadratic equation with respect to B . 
Its solution is as follows 



J3 = A 2 = 



Re±7Re 2 -Re 



2tt 2 S, 



(A.7) 



withRe>Re 



lh ' 



Re, 



is 



from (6), when f5 x 2 = B for 



Re=Re, A , and 



^ 1 =<(p 2 + l)Vl^4^-(l-j ? 2 )^) j^^C 



Sh 



p 2 P 2l P l2 



by =P 22 -pP n ,S= * " " ctgx(p + -) 

(l-p) 71 6, 



S I ^P U P 21 
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Fig. 1 . Family of the curves of the neutral stability (with A, 1 = ) according to (6) and (A.5). In 

Fig. la), a fragment from Fig. lb) is shown. In Fig. la), the upper curve (with two branches) 
corresponds to B =0.463 (for p= 1.527), and the lower one to B =1.099 (for/?=2.239) according 
to (A.5). Fig. la) is overlapped with a figure from [25] (see Fig. 12 in [25]) formally using l/2p= 

aS . In [25], a is the disturbance wave number, and £*is the thickness of the boundary layer 
shift when streamlining a thin plate. In Fig. la), dots and dashed lines correspond to the 
experiment [25], and solid lines correspond to the theory of Schlichting (lower one with the 
Roman digit I), and Lin's (the upper one with the Roman digit II). 

5. Table of values of Re ., and /? obtained by the formulae (A.3) and (6) for p corresponding to 
the local minima of Re .. in the approximate formula (6) 



p 


0o (from 
(A.3)) 


Re th ( from 
(A.3)) 


0o ( from 
(A.5)) 


Re th ( from 
(6)) 


1,527 


0,471 


448,455 


0,463 


442,278 


0,674 


1,124 


680,307 


1,101 


678,482 


0,447 


1,368 


1095,455 


1,358 


1093,824 


2,239 


1,100 


1983,171 


1,099 


1981,838 


2,791 


0,220 


13095,398 


0,219 


13095,285 


0,359 


1,114 


23816,499 


1,114 


23816,488 



.2 I- 



2 <r 



^-f ,. --"" 



A "^ 



Re 



0,807 
0,807- 
0,806- 
0,808 ■ 
0,805- 
0,805- 
0,804- 
0,804- 
0,803 
0,803- 
0,802- 
0,802- 
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Fig. 2, a) Results of the experiments from [22] ] (blue triangles are used for the turbulent spot 
rear front velocity, and white triangles reflect velocity of the leading front) and [24] (squares and 
rounds for the velocity of the rear front). Phase velocity of the wave solutions [8] (numbers 1-5 
denote a degree of the azimuthal symmetry of the running wave). The velocities are scaled by the 
mean flow velocity in the pipe. Phase velocity (A.6) (upper "straight line" corresponds to the 
plus sign in (A.6), and the lower one corresponds to minus) for p~l.53, p =0,471 that 
corresponds to the absolute minimum Re ( i =448.5 according to (A.3). b) Enlarged 

representation of the lower "straight line" from Fig. 2, a) according to (A.6) (the upper branch in 
Fig. 2, b) corresponds to the plus sign in (A.7), and the lower one to minus. 
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